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Abstract: We review recent developments in the calculation of deep-inelastic struc- 
ture functions to next-to-next-to leading order in perturbative QCD. We discuss the 
impact of these corrections on the determination of the strong coupling a s and the 
parton distributions. 
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Structure functions in inclusive deep-inelastic scattering offer the possibility 
for extremely precise determinations of the strong coupling a s and the parton 
distribution functions. The high statistical accuracy of the present experimental 
measurements, and the data expected from the electron-proton collider HERA after 
the luminosity upgrade, demand analyses in perturbative QCD beyond the standard 
next-to leading order (NLO) corrections. In order to match the experimental preci- 
sion, it is therefore necessary to calculate higher order perturbative QCD corrections 
for the structure functions F 2 , F 3 and F L , in particular the complete next-to- next-to 
leading order (NNLO) corrections. This information is not fully available yet. Some 
time ago, the two- loop coefficient functions of F 2 , F 3 and Fl have been calculated Jj]], 
and more recently, they have been completely checked 0. However, for the three- 
loop anomalous dimensions 7^, only partial results are available thus far. These 
include a finite number of fixed Mellin moments [Q, f|, ||, 0], both for F2 and F 3 , the 
large n/-limit [0] of 7^ and 7^, in the latter case only the coefficient of the colour 
factor ti^Ca, and several terms relevant in the small-x limit ||. 



1. Prospects 

To explore the prospects of a complete NNLO analysis in this situation, it has been 
a fruitful strategy to combine all theoretical information presently available together 
with reasonable assumptions about the parton distributions or the functional form 
of the three-loop splitting functions (x) . These investigations resulted in approx- 
imate expressions for PW(x) with negligible residual uncertainties for x>10~ 2 . 
Subsequent studies of the NNLO evolution of parton distributions revealed a much 
reduced scale dependence in comparison to standard NLO analyses. Together with 
recent developments to account for correlated experimental errors in analyses 



and a rigorous statistical treament |TT[ of parton distributions, these efforts seem 
to allow for quantitative estimates of the parton distribution uncertainties. Let us 
emphasize, that this task is of basic importance for precise luminosity measurements 
via W- and Z-boson production in upcoming experiments at hadron colliders, and 
therefore essential for future searches of the Higgs particle or effects of new physics. 

The precision determination of the running coupling a s based on the available 
partial information about NNLO QCD corrections has also been performed. The 
CCFR data for xF 3 from z/iV-scattering |T2| has been analyzed [JT3] using the fixed 
Mellin moments ||, |, |, § for F 3 . The SLAC and HERA data from eP-scattering jU 



have been averaged with Bernstein polynomials to extract a s in a direct fit of Mellin 
moments to experimental data |15f| . The analyses agree within their errors and, 
more importantly, indicate that an absolute error for the strong coupling Aa s <1% 
is possible. In addition, the effect of higher twist contributions, i.e. terms suppressed 



as 1/Q , have been studied |T3| and the NNLO evolution of parton distributions has 



even been used to put bounds on squark and gluino masses [[HJ. Finally, it is worth 
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mentioning that the knowledge of a number of fixed Mellin moments of the three-loop 
coefficient functions |3|, f|, ||, |[ enables analyses even beyond NNLO and allows for 
estimates of the effect of the next-to-NNLO corrections. 

In summary, with the complete NNLO QCD analyses for the structure functions 
F 2 , F3 and Fl a new level of precision is reached in comparing theoretical predictions 
to experimental data. Moreover, the chance arises to address new questions which 
could not be studied before. 



2. Progress 

Progress towards the calculation of the three-loop anomalous dimensions and the 
coefficient functions crucially relies on the ability to perform all necessary loop in- 
tegrations. A promising approach |16], I7J which allows to calculate the integrals in 



Mellin space analytically as a general function of N, is based on the optical theorem 
and the operator product expansion (OPE). The parameters of the OPE are directly 
related to the Mellin moments of the structure functions. For F 2 we can write 

1 

F 2 N (Q 2 ) = J dxx N - 2 F 2 (x,Q 2 ) (2.1) 


= E C»(%,a s )Ai iN (^) , 

i=Q,q,g V V ) 

and similar relations define F® and Fff . Here, C 2 j denote the coefficient functions 
and Ap N the spin averaged hadronic matrix elements of singlet operators O q , s 
and non-singlet operators O a , a = 1, 2, ... , {n 2 — 1), of leading twist. The coefficient 
functions and the renormalized operator matrix elements in eq.( [2.1| ) both satisfy 
renormalization group equations. Due to current conservation they are governed by 
the same anomalous dimensions 7^, which determine the scale evolution of deep- 
inelastic structure functions, 



k=a,q,g ( 



d 




+ ljk (a s (fi 2 ))]A^ N (fi 2 ) =0, (2.2) 



k=a,q,g K 



^W 2+Pia ^ 2)) d^}^ 



V2 



l jk (a s (^))\C^ k (^,a s Gu 2 ) 1=0, (2.3) 



where j = a, q, g and j3 represents the /^-function of QCD. 

The task is the calculation of the coefficient functions C 2 j and the anomalous 
dimensions 7^ in dimensionally regulated perturbation theory. In practice, at a 



2 



given order in perturbation theory, this amounts to calculating the iV-th moment of 
all contributing four-point diagrams with external partons of momentum p,p 2 = 
and photons of momentum q, which is precisely the coefficient of (p ■ q/q 2 ) N ■ 

In order to do so, it is useful to set up a hierarchy among all diagrams depending 
on the number of p-dependent propagators. We define basic building blocks (BBB) 
as diagrams in which the parton momentum p flows only through a single line in 
the diagram. Composite building blocks (CBB) are all diagrams with more than one 
p-dependent propagator. At the three-loop level, there are 10 BBB's, and 32 CBB's 
respectivley, which correspond to genuine three-loop topologies of the ladder, benz 
or non-planar type. In addition, there exist numerous BBB's and CBB's at three 
loops with an effective two-loop topology and a self-energy insertion in one line. 

For the BBB's, the single p-dependent propagator, say l/(p — I) • (p — I) in a 
loop with momentum I, can be expanded into a geometrical sum using p 2 = 0. Then, 
scaling arguments require the final answer for the iV-th moment to be proportional to 
the coefficient (2p-l/l 2 ) N . Thus, one is left with two-point functions with symbolic 
powers of scalar products in the numerator and denominator. For these objects, one 
sets up a reduction scheme, that relates the BBB under consideration to simpler dia- 
grams, where certain lines are eliminated, such that the topology simplifies. For the 
CBB, a straightforward expansion of the p-dependent propagators leads to multiple 
nested sums, which in general are very difficult to evaluate. Hence, one has to seek a 
reduction scheme, that maps a given CBB onto BBB's. This is achieved, if one can 
remove a p-dependent propagator. If one can get rid of a p-independent propagator, 
usually the topology simplifies. 

For the BBB's and the CBB's alike the reduction schemes are determined with 



the help of integration-by-parts identities |18|, [19] and scaling identities [g]. The 
reduction identities often involve explicitly the parameter N of the Mellin moment 
and sometimes one has to set up difference equations in N for the iV-th moment 
F(N) of a diagram , 

a (N)F(N) + ai (N)F(N - 1) + (2.4) 
. . . + a n (N)F(N - n) + G(N) =0, 

where G(N) denotes the iV-th Mellin moment of simpler diagrams. In the reduction 
schemes for three-loop diagrams, we have encountered difference equations up to the 
third order. First order difference equations can be solved at the cost of one sum 
over T-functions in dimensional regularization, D = 4 — 2e. The T-functions can be 
expanded in e and the sum can be solved to any order in e in terms of harmonic 
sums [|(], ^TJ. Higher order difference equations could be solved constructively. In 
general, the approach to calculate Mellin moments of structure functions relies on 
particular mathematical concepts [p2]| , such as harmonic sums p0| , pTf and related 
sums ||23[| . Subsequently, the inverse Mellin transformation to x-space requires har- 
monic polylogarithms 0], or more general, multiple polylogarithms Difference 
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equations for the evaluation of Feynman diagrams, although not in the context of 
Mellin moments, have recently also been studied by other authors |26|, |27 . 

Let us present an example of a diagram of the non-planar type, that gives rise 
to a third order difference equation, 

/ d ° h d% d% i\{ P + h)ni...ii ' (2 ' 5) 

where h = h - q, h = h - h + h ~ q, k = h - q, h = h - h and k = h - h, see 



ref. 1 28] for the conventions of the momentum flow. The result for the iV-th Mellin 




moment of this diagram is given by the coefficient c^, 

= <*(^)"- (2.6) 

The solution for cjy is expressed in terms of harmonic sums of weight six, as it is 
expected for the finite terms of a three- loop diagram. We obtain for cn, 

cn = (2.7) 

(_l)N 

^^-(85-3,-2(^+1) + 85 2 (iV+l)C3 
+4V 3 (iV+l) -4S 2)3 (iV+l) 
+4S 3 ,- 2 (iV+l) +4S 3 , 2 (iV+i) + 10C 5 ) 

+^r( 45 -3,-2(A r +l) +4S_ 3 , 2 (iV+l) 
+8S_ 2 (iV+l)C 3 + 4S- 2 ,_ 3 (iV+l) 
-4S_ 2 , 3 (iV+l) + 85 3i _ 2 (iV+l) + 10C 5 ) • 

This illustrates nicely that our method will not only provide the anomalous dimen- 
sions, which are proportional to the single pole in e in dimensional regularization, 
but also the coefficient functions which are determined by the finite terms at three 
loops. 

In a systematic study we could set up, solve and program the complete reduction 
identities for all three-loop BBB's and all genuine three-loop CBB's of the ladder, 
benz and non-planar type. We have used FORM for this task and performed 
checks at all stages of the calculation with the standard MINCER routine PS|| by 
evaluating the expressions for a number of fixed values of the Mellin moment N. 
Optimization of the program requires the tabulation of several thousand two- and 
three-loop integrals in order to evaluate a given Feynman diagram in reasonable time. 
Thus far, the creation of these tables has already used a large amount of computer 
time. The complete database of Feynman diagrams for the structure functions F 2 , F 3 
and Fl has been generated with QGRAF |Kj and contains roughly 11000 diagrams 
up to three loops. At present, the implementation of the reduction scheme is not 
complete yet. We expect that finishing the program along with the testing and finally 
the actual calculation of all Feynman diagrams still requires a lot of work. 
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3. Conclusions 

The complete NNLO perturbative QCD corrections for deep-inelastic structure func- 
tions significantly reduce the theoretical uncertainties in the determination of the 
strong coupling a s and the parton distributions. A precise knowledge of these quan- 
tities, including a quantitative error estimate, will be particularly important for new 
hadron collider experiments. 

The present approach based on the OPE, to calculate the Mellin moments of the 
three-loop structure functions seems to allow a successful completion. The approach 
relies on the ability to solve all nested sums as functions of N in terms of harmonic 
sums, to set up and solve difference equations in N and, finally, to reconstruct the 
complete analytical expressions of the results in x-space by means of an inverse Mellin 
transformation. 

As far as the theoretical developments are concerned, we believe that there is a 
realistic chance for very high precision measurements in deep-inelastic scattering. 
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